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Abstract 

At nonzero temperature in QCD, about the deconfining phase transition there is a "semi" quark 
gluon plasma (semi-QGP), where the expectation value of the (renormalized) Polyakov loop is less 
than one. This can be modeled by a semiclassical expansion about a constant field for the vector 
potential, Aq^ which is diagonal in color. We compute the shear viscosity in the semi-QGP by using 
the Boltzmann equation in the presence of this background field. To leading, logarithmic order 
in weak coupling, the dominant diagrams are given by the usual scattering processes of 2 — > 2 
particles. For simplicity wc also assume that both the number of colors and flavors are large. 
Near the critical temperature Tc, where the expectation value of the Polyakov loop is small, the 
overall density of colored fields decreases according to their color representation, with the density 
of quarks vanishes linearly with the loop, and that of gluons, quadratically. This decrease in the 
overall density dominates changes in the transport cross section. As a result, relative to that in 
the perturbative QGP, near Tc the shear viscosity in the semi-QGP is suppressed by two powers 
of the Polyakov loop. In a semiclassical expansion, the suppression of colored fields depends only 
upon which color representation they lie in, and not upon their mass. That light and heavy quarks 
are suppressed in a common manner may help to explain the behavior of charm quarks at RHIC. 
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I. INTRODUCTION 



The collisions of heavy ions at ultrarelativistic energies have exhibited a multitude of 
surprising results [1-4]. Experiments at the Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven have shown that in peripheral collisions, the anisotropy with respect to the 
reaction plane, or elliptic flow, can described by nearly ideal hydrodynamics, with a small 
value for the shear viscosity. In kinetic theory, the shear viscosity is inversely proportional 
to the scattering amplitude [5], so a small shear viscosity could be due to strong coupling. 
This suggests that is fruitful to consider the analogy to A/" = 4 supersymmetric gauge 
theories, where one can compute at infinite coupling for an infinite number of colors. For 
A/" = 4 supersymmetry, the ratio of the shear viscosity to the entropy density is l/47r, 
which is conjectured to be a universal lower bound [6, 7]. The shear viscosity has also been 
computed from numerical simulations on the lattice [8], in the hadronic phase [9], and in 
the perturbative quark gluon plasma (QGP) [10-19]. 

In this paper we adopt an alternate approach from A/" = 4 supersymmetry. Instead of 
trying to work down from infinite coupling, we work up from small coupling. While per- 
turbation theory at nonzero temperature is badly behaved, resummed perturbation theory 
works down to much lower temperatures, to a few times Tc, where Tc is the critical temper- 
ature for deconfinement [20-24]. This suggests that nonperturbative effects dominate the 
region near T^, which we have termed the semi-QGP [25-29]. In this paper, the fourth in a 
series [30-32], we compute the shear viscosity in a simple approximation for the semi-QGP. 

The basis of this approach are measurements of the order parameters for deconfinement. 
Consider a straight Wilson line in the direction of imaginary time. 



where P denotes path ordering, T is the temperature, r the imaginary time, and Ar is the 
timelike component of the gauge field, in the fundamental representation. The Wilson line is 
a matrix in color space, and so is not gauge invariant, but its eigenvalues are. The simplest 
measure of the eigenvalues of the Wilson line is its trace, which is the Polyakov loop, 

«=^trL, (2) 

This Polyakov loop is directly related to the propagator of an infinitely heavy quark, which 
for Eq. (2) is in the fundamental representation. Heuristically, one can view the Polyakov 
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loop as measuring the excess free energy fq which arises from adding a colored, heavy quark 
to a thermal bath, {i) ~ exp(— /g/T). 

This Polyakov loop is a bare quantity, but a renormalized quantity can be extracted 
from simulations on the lattice [32-36]. These indicate that the expectation value of the 
Polyakov loop in Eq. (2) is approximately constant and near one for temperatures greater 
than ~ 3 Tc. Below ~ 3 Tc, however, the Polyakov loop is significantly less than one. In 
a pure gauge theory, the Polyakov loop is a strict order parameter for deconfinement, so 
below Tc the excess free energy fg is infinite, and the expectation value of the Polyakov loop 
vanishes identically. With dynamical quarks, the Polyakov loop is only an approximate order 
parameter: because of screening by quark antiquark pairs, below Tc the excess free energy 
fq is finite, and {£) is nonzero. In a gauge theory with three colors and three light flavors, 
as in QCD, numerical simulations find that the Polyakov loop is small at Tc, (£) ~ 0.2, and 
essentially vanishes below ^ 0.8Tc; see, e.g., Fig. (13) of Bazavov et al, [34, 35]. We dub 
the region in which the Polyakov loop deviates from one as the "semi"-QGP [30-32]. 

Similar considerations have also motivated what is known at the Polyakov~Nambu-Jona- 
Lasino model [29] . In these models, the only dynamical variable related to deconfinement is 
the simplest Polyakov loop, that in the fundamental representation, Eq. (2). For an SU{Nc) 
gauge theory, however, the thermal Wilson line has A^^c — 1 independent eigenvalues. These 
are characterized by traces of higher powers of L, trL^ through trL^'=~^. These higher 
moments of the thermal Wilson line are not directly accessible on the lattice, since loops in 
different representations have distinct renormalization constants. These loops have a simple 
physical interpretation, as the propagator of a heavy quark in that representation. 

In the end, we parametrize our results in terms of the Polyakov loop in the fundamental 
representation. Even so, in intermediate steps, we could not represent the computations by 
a simplified model in which only the simplest Polyakov loop enters: the matrix structure is 
essential. We discuss this further in the Conclusions, Sec. V. 

Our approach to the semi-QGP is the following: To the classical Lagrangian, we add 
terms which drive confinement, such as [26] 

■^eflf = &fuzzy T'^ |tr ; (3) 

there are many other possible terms [25-27]. To date, a simple form for the effective La- 
grangian for the semi-QGP has not been obtained; this would allow one to compute both 
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the pressure and the renormahzed loop from the same effective Lagrangian [37]. Clearly 
analysis from numerical simulations, especially in effective models, is essential to gaining 
this understanding [28]. 

We shall see that knowing the full form for the effective Lagrangian is irrelevant to the 
question we address in this paper: near Tc, where the expectation value of the Polyakov 
loop is small, is the shear viscosity suppressed, or enhanced? Naive expectation, based upon 
kinetic theory, indicates that as the cross section decreases, the shear viscosity increases. In 
contrast, we find that the shear viscosity in the semi-QGP is smaller than in the perturbative 



We note that our approach to the semi-QGP breaks strong resemblance to "double-trace" 
deformations of the vacuum theory [38]. In these models, a term such as Eq. (3) is added 
to drive the theory into a confining phase. 

This paper is organized as follows: In Sec. II, we consider a semiclassical expansion of 
the semi-QGP. For technical reasons we assume that both and Nf are large. In Sec. Ill, 
we discuss how to compute transport coefficients in the semi-QGP using kinetic theory. 
This is a standard approach, except that kinetic theory is in the presence of a background 
field, Ar in Eq. (2). This background field is important near the phase transition, and sup- 
presses the shear viscosity. Section IV gives numerical results. In Sec. V we discuss possible 
phenomenological implications of our results. We find that in a semiclassical expansion of 
the semi-QGP, near all colored particles are suppressed in a universal manner. This is 
natural, as a type of "bleaching" of color as T — t- T+. In particular, this color bleaching 
is independent of the mass of the field. We suggest that this might help to explain the 
otherwise puzzling results on the behavior of charm quarks at RHIC. 

II. EFFECTIVE THEORY IN THE SEMI-QGP 

A semi-QGP is characterized by the Polyakov loop, Eq. (2). Static quantities, such as 
the pressure, are determined by an effective Lagrangian of Polyakov loops [25-29]. What 
we require is an effective theory in real time, which can be used to compute transport 
coefficients. Start with the partition function for QCD, 



QGP. 




(4) 
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where £ is the Lagrangian, 

C = -]^tiFl + ij{ip-m)^, (5) 

with F^^ = dfj,Ay — dyAfj_ — ig[Afj_, Ay]] ip and ip is the quark field, If) = 7^(9^ — igA^), and 
m is the quark mass. C denotes a path in complex time, Fig. (1). We work in Minkowski 
spacetime with a metric, g^y = g^^ = diag(l, —1, —1, —1). The timelike gauge field on the 
imaginary time in Eq. (1), A^, corresponds to iAq in this notation. We defer a discussion 
of a fully self-consistent approach, where terms such as Eq. (3) are added, to later analysis. 
Including such a term would not alter our results at leading order in g, because to leading 
order the shear viscosity is dominated by the scattering of hard particles, with momenta 
~ T, while terms at Eq. (3) affect fields at soft momenta. 

For the generators of SU{Nc), we use what may be called 't Hooft's double line basis [31, 
39]: 

[t%<i = ^kSd-^S'^'S^<i)- (6) 



Somewhat unconventionally, we denote indices in the fundamental representation as 
a,b,c,d, . . . = 1 . . . Nc- In the double fine basis, generators are denoted by a pair of in- 
dices in the fundamental representation, ab, so that in all there are generators. This is 
one more generator than in an orthonormal basis. As an overcomplete basis, the trace of 
two generators is not a delta function, but a projection operator: 

1 1 

-j-j, ^ab ^cd _ pab,cd pab,cd ^ad ^bc ^ab ^cd 

The commutator of two generators is 

= I ^ fiab,cd,ef)^fe^ ^g^ 

where f'y°-^'"'-^'^f) is the structure constant, 

j{ab,cd,ef) _ _J_ ^^ad ^cf ^eb _ ^af ^cb ^ed^^ ^ ^g-^ 

The double line basis is obviously convenient at large Nc, as then the terms ~ 1/Nc in 
Eqs. (6) and (7) can be dropped. It is also of use at finite Nc, however, since then each line 
in a generator represents a flow of color charge, and charge conservation is obvious; this is 
especially true in the presence of a background field. [31]. 



We denote four momenta as P'^ — {p^,p), where po is the timehke component of the 
momenta, and p the spatial momentum. In thermal equilibrium at a temperature T, po — 
iujn, where a;„ is a Matsubara frequency: a;„ — 2mrT for bosons, and = {2n + l)7rT for 
fermions, where n is an integer. In imaginary time, a nonzero value of the Polyakov loop is 
modeled by taking a constant field for Aq, 

[AXb = iSab y . (10) 
We expand about the classical field, A^, in fluctuations, B^: 

= V^O (11) 

The covariant derivatives in this particular background field are especially simple: 

[i{do - igA'Ma {po + iQ'')A ; 

[t{do-tg[A''^,.])B,u^iPo + tQ''')[B^U ; g'^'=g'^-g^ (12) 

Note that the adjoint covariant derivative involves the quantity Q"'', which is the difference 
of two color charges. The corresponding statistical distribution functions are 



where I • • • ji^ is defined as 



, , , +2; Re 2; > 

l^k = { ■ (14) 

-z Re z <0 



From the statistical distribution functions we see that the background field acts like an 
imaginary chemical potential for color. As such, in a given field these distribution 
functions are complex, and so unphysical. The only physically meaningful quantities are 
integrals over distributions of Q, and these give results which are sensible. 

To see how this comes about, consider summing a quark propagator over its color indices. 
This enters, for example, in the computation of the pressure at leading order. The sum is 

]J;E^<5=«WT-E(-r'e-'^^trL". (15) 

a n=l 

We can always assume that the vacuum expectation value of tr L is real. In the pure glue 
theory, this can be enforced by a global Z{Nc) rotation; with dynamical quarks, this is 
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FIG. 1: Complex-time path. 



automatic. In Eq. (15) there is an infinite series in powers of trL", the expectation values 
of which are real. The computation of the shear viscosity in this paper provides another 
example where integrals over Q's give physically sensible results. 

We wish to compute correlation functions near thermal equilibrium, and so need to 
continue from imaginary to real time. We adopt the usual path, illustrated in Fig. (1), 
along C = Ci U C2. This includes integration in real time. Ret, along Ci from an initial 
time ti, to a final time, t/, and back (in practice, both times are assumed to be infinite). 
Then one integrates in imaginary time, along C2 from Imt : — i- — 1/T = — /3. 

It is natural to take the background field only for the part of the path in imaginary time. 
Consider an ordinary chemical potential, introduced as a Lagrange multiplier for the number 
operator. The chemical potential does not affect the Hamiltonian, nor the evolution in real 
time: it enters only to change the statistical distribution functions in thermal equilibrium. 
Thus we do the same for an imaginary chemical potential for color: along Fig. (1) we take 
the background field to be nonzero along C2, and to vanish along Ci. 

While we will not use the real time formalism, it helps to understand the choice of 
background field in time. The propagator is a two by two matrix. 



ab,cd 



dfiu Pab,cd Dab{K) 



(16) 



where 




(17) 
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with 



D^^{K) ={nab{ko) + e{-ko)) 271 6{K^) , 
DlliK)={natiko) + eiko))27r 5{K^) , 



(18) 



in Feynman gauge [40], where 6{k^) is the Heaviside step function. This form shows that 
the zero temperature parts of the propagator, 1/ {K^ ± ie), are independent of the Q's. The 
color charges only enter through the statistical distribution functions, nab{ko). 

Instead of using the real time formalism, we compute in imaginary time, and then ana- 
lytically continue. The color dependent frequency Un becomes 



where ± corresponds to retarded and advanced propagation. 

After computing diagrams, we need to sum over distributions of the Q's. To avoid having 
to deal with trace terms at finite Nc, we work in the limit of large Nc. The pressure is of 
order ~ A'"^ in the hadronic phase, and ~ A^^ in the deconfined phase. Thus a mean field 
approximation, as in Eq. (10), should be a reasonable approximation in the deconfined phase 
at large Nc- At large N^, color sums can be replaced by integrals: 



where we introduce the eigenvalue density p{(f) = da/d(p, with (f = Q/T. Moments of the 
Polyakov loops are Fourier transforms of the eigenvalue density, 

1 r , . . 



In particular, we write the the Polyakov loop in the fundamental representation, £i, simply 
as i. In the confined phase, all moments of the Polyakov loop vanish, in = 0, and the 
eigenvalue density is constant. In the perturbative QGP, the moments of all Polyakov loops 
are close to one, = 1, and the eigenvalue density is close to a delta function. 

Since we do not now have a complete effective theory of the semi-QGP, we do not know 
the exact eigenvalue density as a function of temperature. We do know the expectation 



(19) 




(20) 




(21) 
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value of the simplest Polyakov loop, £, as a function of temperature, and so we consider two 
representative forms for the eigenvalue density. One ansatz is to take a step function, 

Pstep(</',A) = ^^(A-|</.|), (22) 

where A is a parameter characterizing the eigenvalue density. With this ansatz, the expec- 
tation value of the Polyakov loops are 

(23) 

The other ansatz we take from a soluble model in two dimensions, that of Gross and 
Witten [41]: 

, (1 + Acoso9) A < 1 



cos 



2 



A>1 



where \(p\ < 2sin^^ (1/A)^''^ is satisfied for A > 1. For A < 1, £±i = 1/(2A), while all others 
vanish. With this ansatz, there is no simple expression for when A > 1. 

Rather unexpectedly, we shall see that our results for the shear viscosity are almost 
independent of the choice of these eigenvalue densities. We do not know if this is a property 
special to the shear viscosity, or if it is generic. 



III. TRANSPORT COEFFICIENTS 



In this section we consider transport coefficients in a semi-QGP. Hydrodynamics is de- 
scribed by evolution in time of locally conserved currents, such as the energy-momentum 
tensor, which satisfy d^T^^ = 0. For simplicity, we assume that the only conserved charges 
are those for energy and momentum. In hydrodynamics the energy-momentum tensor T^" 
can be expanded as [42] 

T^"" = Su^u" + VA^"" + n^"" , (25) 

where u'^ is the local velocity, u^m'^ = 1; S = T^'^u^Uy is the energy density and V is the 
local pressure, both of which are Lorentz scalars; lastly, A^'^ is a spatial projection operator, 

A^''' = u'^u'' - g^"" . (26) 

For the metric we take g^u = g^^ = diag(l, —1, —1, —1). 
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jjfj^f measures the deviation from thermal equihbrium. Near thermal equilibrium it can 
be expanded in a gradient expansion, 

n'"' = na>"' + C Z\^M"'^ dc^up + ■■■ , (27) 

where 

^ ^,p^ua (^Q^^^ ^ Q^^^ _ 2 A'^f'd^u^ . (28) 

The first order transport coefficients C and ri are, respectively, the bulk and shear viscosities. 
The transport coefficients reflect microscopic properties of the theory, and are an input into 
a hydrodynamical analysis. 

These transport coefficients can be computed from the Kubo formula, which relates them 
to the correlation functions [43]: 

=1 hm - / rfW-*^(t)(K,(x), 7r^^(0)]) , (29) 

C =\ lim - / ciW-*^(t)([P(a;),P(0)]) , (30) 

where V{x) = —T\{x)/3 and Hijix) = Tij{x) + gijV{x). These correlation functions are not 
simple to compute in perturbation theory, though. They are sensitive to pinch singularities, 
and it is necessary to resum an infinite series of diagrams. 

Another way of computing these transport coefficients is to use kinetic theory and the 
Boltzmann equation. To leading order, this is equivalent to a resummation of the relevant 
diagrams in the Kubo formula [13, 14]. We will use the Boltzmann equation to compute the 
shear viscosity in a semi-QGP in Sec. IIIB. Our analysis is a straightforward extension of 
the computation of Arnold, Moore, and Yaffe [16], of the shear viscosity in the perturbative 
QGP, to the semi-QGP. 



A. Shear viscosity in the semi-QGP 

Before going to the technical details of the computation, we give a qualitative discussion 
of how the shear viscosity arises in kinetic theory in the QGP, and how it differs in the 
semi-QGP. 

For the transport theory of a relativistic gas [44] , the shear viscosity rj is given as 

4 

V^ — npX^ip, (31) 
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where n is the number density, p is the mean momentum, and Amfp is the mean free path. 
For hght particles, where the mass is much less than the temperature, the mean momentum 
p ~ T, and the number density is n ~ T^. The mean free path is A^fp ~ 1/ (na), where a is 
the transport cross section, so that the shear viscosity is 77 ~ T/a. 

In QCD, the transport cross section is cr ~ g^ln{T/mY)ehye)/T'^, where mDebye ^ 
\fghijT ~ gT is the Debye mass. There is a Coulomb logarithm, ln(T/mDebye), due an 
infrared singularity from forward scattering [45], so in all. 

The Boltzmann equation gives the same powers of g and In (1/(7), but a strictly perturbative 
analysis gives the wrong overall coefficient. To get the correct result perturbatively requires 
the resummation of an infinite set of diagrams [11, 16]. 

In weak coupling, a is small and the shear viscosity is large. This sounds counter intuitive, 
but is not: when perturbed, it takes a long time for a weakly interacting gas to move back 
towards thermal equiblibrium. 

Conversely, the most natural way to obtain a small shear viscosity is if the cross section, 
and so the coupling constant, are large. The best illustration of this is a A/" = 4 super- 
symmetric gauge theory, where in the limit of infinite coupling and an infinite number of 
colors, the ratio of the shear viscosity to the entropy density, s, is r^/s = l/47r [6]. This is 
conjectured to be a universal lower bound on rj/s. 

The analysis in a semi-QGP is distinct from ordinary kinetic theory. It is characterized 
by the partial ionization of color [30-32], so that both the density of colored particles, and 
the interaction cross section, depend nontrivially upon the Polyakov loop. 

Consider first the density. From Eq. (13), a background Q field acts like an imaginary 
chemical potential for color. For ordinary statistical distribution functions, the only way to 
obtain a small density of particles is if they are heavy, and so Boltzmann suppressed. A 
nonzero Q field provides another. When Q 0, the statistical distribution function for a 
given particle is complex valued. Thus while the distribution functions are of order one in 
magnitude, because of cancelling phases, they can vanish after averaging over the Q's. This 
happens in the confined phase at Nc = 00, where the expectation value of trL" vanishes for 
all n. For a quark, which lies in the fundamental representation, the trace of the statistical 
distribution function is given in Eq. (15), and vanishes term by term. Just above Tc, where 
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i is nonzero but small, the distribution function for such a field vanishes as a single power 
of the loop, 

(^n,)~iV,T3£; (33) 

a 

the bracket (■ ■ ■ ) denotes the average over the Q distribution, and an integral over the 
particles three momentum. 

Gluons in the adjoint representation carry two Q charges, so the Boltzmann expansion 
of the statistical distribution function involves exp(i((5'^ — Q^))- Summing over the color 
indices, Xla^**^" X^b^"*'^'' ~ ^^-^ tiL\ so the density of gluons then vanishes as the square of 
the loop, 

(£n^,) - Nl e ■ (34) 

ab 

we assume that the expectation of ^ is real. This suppression of the density has no analogy 
in the perturbative regime, where the density of massless fields in necessarily a pure number 
times T^. 

The average mean free path is the ratio of the density, divided by a transport cross 
section. For the scattering of two gluons, 

Ai^fp - ^^^'""'"^ r ; (35) 

\l^ab,cd ^o.b ricd (Tab;cd) 

o'abicd is the transport cross section for 2 — )• 2 scattering between gluons with color indices 
ab and cd. 

For the sake of argument, assume first that there are no correlations between different 
colors. If true, then the cross section would factorize, 

( nab ricd CTab-cd) = rLab)"^ {crab;cd) ~ CT t . (36) 

ab,cd ab 

Here cr is a typical perturbative cross section, ~ g'^\og{l/g). The mean free path would then 
diverge at small i, A^^p ~ l/(o"£^). This factor of would cancel a similar factor in the 
density, so that the shear viscosity would be as in ordinary perturbation theory, 77 ~ T/o", 
(approximately) independent of i at small i. 

As we shall see in the next section, however, there are strong correlations between different 
colors, so that the cross section does not factorize as in Eq. (36). Instead, it vanishes only 
quadratically with the loop, 

nab ncd crab;cd) ~ <y , (37) 



ab,cd 
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At first sight this dependence on ^ seems counter intuitive: the scattering process involves the 
scattering of two gluons, with color indices ah and cd. Expanding the statistical distribution 
function for one gluon brings in powers of exp(i((5'^ — Q^))-, while expanding that for the 
other gluon brings in powers of exp(i((5^ — Q'^))- Summing over all color indices, the naive 
expectation is that the cross section is proportional to the fourth power of the loop, one 
from each factor of exp(i(5). The scattering processes, however, include a planar diagram in 
which two of the color indices are equal, say b = c. For this term, the factor of exp{—iQ'') 
from one gluon cancels, identically, against the factor of exp{+iQ'^) from the other gluon. 
This explains why Eq. (37) is proportional to i"^, and not to i"^, as in Eq. (36). 

In the pure glue theory, then, from Eq. (35) the mean free path is approximately inde- 
pendent of the expectation value of the loop, A^^p ~ 1. Even in the confined phase, the 
mass dimension of Q is set by the temperature; thus the typical momentum remains ~ T, 
and p ^ T. Consequently, from Eq. (31) the shear viscosity is proportional to the overall 
gluon density, and so vanishes quadratically with the loop, ~ 

With dynamical quarks the counting differs, although somewhat unexpectedly, the result 
is similar. The cross section for the scattering between a quark, with color a, and an 
antiquark, with color b, is 



If we could factorize the cross section, as in Eq. (36), then this would be proportional to 
the density of quarks, squared, or ~ There is such a contribution when the colors are 
unequal, a ^ b. In addition, though, there is a planar diagram from when the colors are 
equal, a = b. In this case the chemical potential for the quark, exp(zQ"), cancels identically 
against that from the antiquark, exp{—iQ''), so that Eq. (38) is not ~ but constant. The 
mean free path for quarks. 



is then proportional to the quark density, ~ i. The average momentum remains ~ T, so 
that the shear viscosity is the product of the quark density, ~ £, times the mean free path. 
Thus like gluons, in all the quark contribution to the shear viscosity is r/'^'^ ~ 

Putting in the factors of the coupling constant, we find that either in the pure glue theory, 
or in the theory with dynamical quarks, that the shear viscosity vanishes quadratically in 




(38) 



a,b 



qk 
mfp 




(39) 



13 



the loop, 



T3 



(40) 



relative to the result in the perturbative QGP. 

We show later that the results for quarks is in fact numerically close to that of the pure 
glue theory. While we consider a large limit in which Nf is also large, this similarity 
between the shear viscosity in the pure glue theory, and with dynamical quarks, is rather 
unexpected. 

B. Derivation of the shear viscosity in the semi-QGP 

In this section, we use the Boltzmann equation to compute the shear viscosity in the 
semi-QGP [10-12, 16-19]. Of course we assume that the coupling constant is small enough 
so that a quasiparticle picture is applicable, and kinetic theory is valid. The computation is 
precisely like that in the ordinary QGP, except that we need to compute in the presence of 
a background field for Aq. In this section we consider a purely gluonic theory, considering 
the case with dynamical quarks in the next section. 

If the statistical distribution function for a gluon with color ah is fab-, then in kinetic 
theory the energy-momentum tensor is 



The factor of two in Eq. (41) comes from the following relativistic normalization of the phase 
space including color and spin: 



where s denotes the spin of the gluon, and E = \p\ is the energy of the gluon. We are 
interested in a system slightly away from local thermal equilibrium, so we write the statistical 
distribution function as 




(41) 




(42) 




(43) 



where /^^ is the gluon distribution function in local thermal equilibrium. 



f(0) 



1 



(44) 
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FIG. 2: Feynman diagram of the gluon scattering with the normal and double line notations in 
the t-channel. 

with pt^'"-^ = pt^ — iS'^^Q"''', and 6 fab measures the deviation from equilibrium. 

The constant background field A'^ induces an imaginary color charge. We assume that 
after averaging over the Q's, that the net color charge vanishes. A spatially dependent color 
charge arises in related problems, such as for 't Hooft loops, when there is a background 
color electric field [31, 46-48]. 

As discussed in Sec. II, the background field vanishes along the part of the contour in 
real time in Fig. (1). Thus the time derivative in the Boltzmann equation does not involve 
the background field, 

2P>^dJaM^,p,t) = -Ca,bAf]- (45) 

We do not consider a force term, which may be dynamically generated in an expanding 
plasma [17]. In Eq. (45), C[f] is the collision term for the scattering of 2 — )■ 2 gluons. 



4 r 

CaMf] =n / rfPi(2vr)V(Pi+P2-P3-P4)|A^r 
i=2 

X (^/ai 61/0262(1 + fasba)^^ + fa4bi) ^ (1 + /ai6i)(l + f a2b2) f azb^ f a^b^ 



(46) 



where is the scattering amplitude, which depends upon the background field. To leading 
order, it is necessary to include not only the scattering of 2 — )• 2 particles, but also processes 
involving splitting, of 1 — t- 2 particles, and joining, for 2 — 1 particles [16, 49, 50]. The 
processes of soft multiple splitting and joining are known as the Landau-Pomeranchuk- 
Migdal (LPM) effect [51]. In this paper we compute to leading order in both the coupling, 
g^, and logarithms oi g. At leading logarithmic order, the result is dominated by the 
scattering of 2 — )■ 2 particles in the t-channel, and the LPM effect can be neglected. 

In the pure glue theory, the corresponding diagram in the t-channel is shown in Fig. (2). 
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In the double line basis, there are four diagrams. The infrared singularity is cutoff by using 
a dressed propagator in the Hard Thermal Loop (HTL) approximation, in the presence of 
the background field [31, 52], 



(47) 



and are the propagators for the longitudinal and transverse gluons in the HTL 
approximation, respectively. For the amplitude at tree level, gauge dependent terms can be 
ignored. Projection operators for the transverse and longitudinal directions are 



fc2 



while the transverse and longitudinal propagators are 



(48) 



i 



2 x-l\ 



ab,cd 



A;2 + 2m2(l- f In^^' 

^ 2 x—1' 



(49) 



J ab,cd 



where x = kQ/\k\. The thermal mass, [m^ 



\ab,cd} 



depends upon the details of the color 



distribution [31]. In general, it is not diagonal, and different components of [m'^]aa,bb, with 
no summation over the indices a and b, mix with one another. In the limit of large N^, 
[nT''^]ab,cd becomes diagonal, and mixing is suppressed by l/Nc- In the following, we work 
only to leading order of Nc, in which the propagators are diagonal: 

DlU^) = Sad he Dl{K) 

= Sad Sbc 
Dab,cdiK) = Sad Sbc D^,,{K) 
= Sad Sbc 



P + 2[m2]„,(l-|lnf±f 



where 



g'NX' 1 

6 2Nr 



(50) 



(51) 



e=l 
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FIG. 3: Kinematics of the two body scattering. 



and 

^(g) = l-6g(l-g), (52) 

where q = Q/{2-nT), and q is defined modulo 1 [31]. 

We note that at large A^c, the combination g'^Nc, and so the Debye mass, are of order 
one. In absence of the background field, Q = 0, ['m'^]ab coincides with the ordinary thermal 
mass, [m\b = 9^N,T^/Q. 

For the shear viscosity, we find that the result is proportional to that in zero field, 

V = r/port ^(^) , (53) 
where ?7pert is the usual result in perturbation theory [16], 

constant of order one, which depends upon Nc and Nj. 
Note that at large N^, the free energy and so the entropy are large, of order N^. The 
shear viscosity is of the same order, so that the ratio of the shear viscosity, to the entropy, 
is of order one. Our task is now to derive TZ{i). 

1. Kinematics and Scattering amplitude 

Before solving the Boltzmann equation we review the kinematics of scattering processes 
for 2 —> 2 particles, involving the exchange of a soft gluon [11]. In terms of the Mandelstam 
variables, 

s = {P^ + P^f , t={P,- P,)' , u={Pi- P,f . (55) 
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The dominant kinematic region in Fig. (2) is |t| ^ |s| ^ as then the amphtude, squared, 
is enhanced by ~ We introduce the quantities: 

K K K K 

in which the Mandelstam variables become 

s = {P + P')\ t = K\ u = {P-P'f. (57) 
There are three angles which enter, as shown in Fig. (3): 

cos 6 = k ■ p , cos 6' = k ■ p' , cos cp = — 1 — — — :r-^ , (58) 

where unit vector are written a.s p = p/\p\, etc. For forward scattering, 6 is close to 6', and 

p ■ p' — 

cos if ~ _ ^ , (59) 
using cos^' ~ X = A;o/|fc|. The Mandelstam variables become 

s ~2pp'(l — x^)(l — cos(f) , 

U — 2pp'{l — X^)(l — COS<y9) , 

t-~k\l-x^), (60) 
with p = \p\, p' = \p'\, and k = \k\. The amplitude in the t-channel is 



M = JZ^iPi^Ps) D,,,,,^,,{K) J:^i{P2, P,) , (61) 

where J'^'^!' is the color current, and D^yab^cd is the gluon propagator in Eq. (47). For forward 
scattering. Pi ^ P3 P, and the form of the current is dictated by Lorentz symmetry, 

jr;^/3(Pi,P3)^2z^P^tf 5,,,3, (62) 

where tf' is a generator of SU{N) in the representation r. For the adjoint representation. 

The spin or helicity is denoted by Si and S3. The amplitude for gluon scattering becomes 



ab,cd 

= -4g' 6s„s,6s,,s, ^^13 ^24 PP \Dt{K) + {1 - x') cos ^ D^\K) 



ab 
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using the shorthand notation tfj = [t^b.jaif)i,ajbi_ ^]^g square of the amphtude is 

ab,cd 

X ^Dl\K) + (1 - x^) cos ip D^{K) Dl^{K) + (1 - x^) cos cp D^^{K) 



(65) 



In bare perturbation theory, the thermal mass vanishes, and this becomes 

\M\'' = Ag'Y. til t\\ tli 41 5,,,3 5,,,,,^ . (66) 

ab,cd 

It is rather useful to write phase space in terms of the variables k, p, p', x and ip, instead 
of pi to p4. By energy-momentum conservation, 

6{Ei + E2-E3- Ei) ^6(k{cose - cos^')) 

1 (67) 



-5(cos6' — cos^^') 
k 



the integral over phase space becomes 



fl I dm2n)H^'\P, + P2-P3- P4) 

1 ^ /■ k 

~ 22(27r)6 n 5Z J dkkdpp dp' dx dp {27rf6^^^ {pi -p- -) . 



(68) 



i=2 aibi,Si 



2. Chapman- Enskog method 

Because the Boltzmann equation is nonlinear, it is not easy to solve exactly. Instead, 
we use the Chapman- Enskog method to estimate the transport coefficients [53]. Start with 
a solution of the Boltzmann equation local thermal equilibrium /^°\ Eq. (44). By the 
conservation of energy-momentum and charge, the collision term vanishes, and gives the 
following relation: 

AO) AO) M) AO) 

J aibi J a2b2 •'12363 ^0464 (69) 



^ ' Jaib2 ^ ' Ja2b2 ^ ^ J a-^b^ ^ ' ^0464 

We are interested in a nonequilibrium system which is close to local thermal equilibrium, 
and so expand the distribution function around fj^l\ 

u-fS + ^fT + ^'fT + ---- (70) 
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This corresponds to expanding in spatial derivatives, or equivalently, the Knudsen number. 
The parameter e keeps track of the order of derivatives, and is taken to one after computa- 
tion. /*^") gives the transport coefficients of n-th order, e.g.^ f^^^ gives the shear and bulk 
viscosities. The collision term is expanded in powers of e as 

oo 

C[/] = ^e"C("), (71) 

n=0 

where 

1 ^" 

(72) 



1 r/" 

^(n) ^ " ^rji 



e=0 

(0) 



The first term, C^^\ obviously vanishes, because /^^ is the solution in thermal equilibrium. 
In the left hand side of the Boltzmann equation, the spatial derivative, (9j, is taken to be of 
order e. In the Chapman-Enskog method, the time derivative is expanded in powers of e, 

oo 

9, = ^6"9f). (73) 

n=l 

We show in the following that Sf"'' is determined by energy-momentum conservation and 
local thermal equilibrium, Eqs. (81) and (82). Z^^"* is parametrized as 

/i? = /r(l + /r)^^"^ (74) 
for useful. To e", the Boltzmann equation is written as 

2P''5j^)/„^,^ = -£^W (75) 
where the partial derivative, to order e", is 

d^^'^fa, =p'j2 ^J'Vir'^ + P' d. ft'^ , (76) 
fc=l 

and 

Kit = -C!:l + C <!>^:l . (77) 
We introduce the linearized operator £, 

Ml = n / (2vr)^5(^)(Pi + P,-P,- P,) \M? 

i=2 

X /il /SI (l + /il) (l + /II) + '^t' - - 4"^ 



(78) 
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In vanishing background field, the operator L is self-adjoint and positive semi definite: 

y dr <P^ iC<P) = j dr{C <py<P > for <P^0. (79) 

While positivity is not automatic in a nonzero background field, we assume it is, after 
averaging over color. Our analysis shows no signs of such unphysical behavior. 

Equation (75) is a linear equation for <p("). Since Kn and d^ff* fab are as functionals of 
for < < n — 1, and does not include fl^\ can be obtained order by order. 

There is an ambiguity, however. C contains zero eigenvalues, which correspond to con- 
served currents, ^ = and 1. To remove this ambiguity, we require that at each order in 
e, the deviation from thermal equilibrium is orthogonal to the zero modes, 

dr P° /f;;^ = (80) 



for n > 1. Then the conserved charge density only depends upon /^||\ 

= 2 j dP P° ff^ . (81) 

We also require that energy and momentum are conserved order by order in e, 

2 j dPP'^P'f^l'^ + 2diJ dPP'P'f^"^-'^ = , (82) 

which defines the time derivative to e"^. 

The viscosities are determined at leading order, 

2P''5i'Va... = -/:<L, (83) 
after using Ki = 0. The left hand side is 

2P^ fab = -2 P^'-^^ P^ fif (l + /I?) {u, + /3 9«n.) . (84) 
Working in the local rest frame, where uq = 1, Ui = 0, dl^\o = and d^^^i ^ 0, 



(85) 



3 

hj 

By assumption, the last term, proportional to Q"^, vanishes after averaging over color, and 
so we drop it in computing the shear viscosity. This equation contains which is yet to 
be determined. From Eq. (82), to leading order 

j dP P'^P^d^^f^S = (86) 
21 



which reads 



P di'Ki + ^^/3 = 0. (87) 

Using Eq. (87) in Eq. (85), only last term on the left hand side survives, 

2 d^^f = -/3 (l + pY a., . (88) 

The trace of this equation vanishes, so that the bulk viscosity vanishes in this calculation. 
At linear order, the Boltzmann equation becomes 

+ J>Vcr,, = £<i^. (89) 

Acting with the inverse of the collision operator, on both sides of Eq. (89), we obtain 
the formal solution. 

Note that the condition of orthogonality, Eq. (81), ensures that zero modes are projected 
out, so that the inverse operator is well-defined. We parametrize ^^^^ as 

^>(i) = (pV - ^S''^ a,, xip) = X''{P) ■ (91) 

Remember that aij is traceless, so that here we can add p'^6^^aij/3 on the left hand side, 
and take x*-' to be traceless. The transport coefficients are independent of one another, so 
we write 

Pfil (l + /il) - ^^V) = Cx'l,^ . (92) 

Using the solution of the Boltzmann equation, we rewrite il*-' in Eq. (25) to obtain 

i7(^)^^- = 2jdr pV /f;) = 2 1 rfr pV /S^ (i + /S^) <p^'^ , (93) 

where il^^)*-? denotes e^. Using Eq. (91), and performing the integral over angles gives 

J ^ pVp¥ = i^''^'' + + ^''^'') ' (94) 

we find 
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The rank tensor two in the integrand is nothing but Cx^^ in Eq. (92), so that 

21 

~5 



IT 



I dr {x'^Cx'') . (97) 



The shear viscosity is then 

2T 

In order to determine Xiji we need to know the inverse of a hnear operator, with an infinite 
number of matrix elements. 

We first find a formal solution, expanding x in orthogonal polynomials, x = Sc„Xn(p), 
and 

j dr fab{l + fab) P'^ Xn{p) Xm{p) = dn 5„m , (98) 

where dn is a normalization constant. We choose Xn = Ylm=o ^nmip/T)'^ to be a monic poly- 
nomial, i.e., the leading coefficient is unity, hnn = 1- This polynomial can be systematically 
obtained by Gram-Schmidt orthogonalization. The first three terms are 

Xo=l, (99) 
XI = (I) +610, (100) 



X2 = [^) +&2l(^j+&20, (101) 



where the bnm are dimensionless constants 

ho 



bio = -^, (102) 



hih-^ — hi 
hoh2 - h\ 



^20 = rf , (103) 



hih2 — h^h 



where 



h^ = l3^-n j dr fabil + fab) P 

_ r{6 + n) g |trL^' 



4+n 



(105) 



27r2 ^ A;5+" 
fc=i 



The normalization constants are obtained from Eq. (98)-(101) as 

do = ho, (106) 

di = h2 + hiho , (107) 

d2 = hi + h^h2i + h2h2Q . (108) 
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Taking J dF x^i to Eq. (92), the Boltzmann equation reduces to a matrix form, 

Sn ^ ^ l~'nm i (109) 

m 

where 

Sn = ^ jdr ^ fil^ (i + fil^) (^pY - ^6^^^ , (110) 

with 

Cn^ = j dr (x^Cx'ri) . (Ill) 

Formally, the solution can be written as 

x'^ = j2xnc''usm. (112) 



Inserting this into Eq. (97), we obtain 

~5 



v = ^-^J2s4^-'L^s^- (113) 



n,Tn 



From Eqs. (94) and (98), 

2 

Sn = do 6n,0 , (114) 

from which the shear viscosity is 

•)=^ ("5) 

This involves the zero- zero component of the inverse of C. In general, finding the inverse of 
C is not easy, as it involves an infinite number of matrix elements. It is thus necessary to 
truncate the number of matrix elements. This gives a lower bound on the solution, because 
Eq. (113) has a quadratic form, and by Eq. (79), £ is positive semidefinite. We will see the 
solution converges quickly as the size of the matrix for C increases, with the solution for 
one matrix element within 0.6% of the exact solution (see Fig. (4) and related discussion 
following Eq. (137)). Using just one matrix element, we obtain an approximate form for the 
shear viscosity. 



We now turn to estimating the matrix element, 

=\tlf dr. (2vr)^5(^) (P, + P. - Ps - P,) \M\^ C. (l + fll) (l + C 

1=1 

(117) 
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The coefficient 1/4 is caused by symmetrizing the colhsion term. We are interested only 
in the forward scattering amphtude, where the exchanged momentum, k, is small. In the 
integrand of Eq. (117), to leading order in k, 

Xn'(pi) - xiiips) = xii (p + ^ ) - xii (p ^ 



2 y 2 

p^k^ + fcV - 2{p ■ k)^^ xnip) (118) 

+ (pY - ^) (P ■ ^ {p'Xnip))' + Oik') . 

Further, 

(xn''(pi) + Xn'(P2) - x^iips) - xliip^)) (xZiPi) + xtiP2) " xZiPs) ' xZiPi)) 

2 1 

~ 2 [p'P - {p ■ kf) Xn{p)Xm{p) + -{p-kf — {p\n{p))' {P^Xm{p))' + {p ^ p) 

6 p 

= k'(2{l - X') p' Xn{p)Xm{p) + ^ {P\n{p))' {p\m{p))') + {p ^ p') 
= k''Fnm{P,P',x) . 

In the second line, we dropped the cross term (XniPi) ^ Xn{P3)){Xm{P2) — XmiP^)) that 
vanishes in the integral of 6 and (p. The function Fnm{p,p' ,x) is independent of k; the 
matrix element becomes 



(119) 



1 ^ r 

J^nm^jll dn (27r)^5(^)(Pi + P2-P3- Pa) \M\' 
i=i ^ 

X fit c. (1 + /SI) (1 + /a) k' fup,p', X) 



(120) 



Expanding the statistical distribution functions at small kQ, 

00 00 00 00 



C, C (1 + /S.) (1 + C) = E E E E -p [ - /^E'vf^. - ''Q' 



rai=l n2=l n3=0 n4=0 j=l 
00 00 00 00 



EEEEew. 



where 



~ > > > > y[LJ,n] 

711=1 n2=l n3=0 n4=0 

X e-^('^i(P+'¥)+"2(P-'?)+"3(p'+1l)+n4(p'-1i)) 



4 



(121) 



g{Q, n) = exp ^ n.Q"^''^] . (122) 
Note ^(0,^) = 1. Equation (120) has no singularity at ko = 0, so taking fco — )■ 0, 



00 00 00 00 



C, (1 + /SI) (1 + /il) - E E E E „)e-««".«'>'+i»-«)'''> . (123) 



ni=l n2=l n3=0 n4=0 
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Writing the collision term using this, and using the explicit form of the squared amplitude 
in Eq. (65) as 

-A 



C =2 



9 



-rt CO CO CO CO „QQ 

n E E E E E ^(^' ^) ^13 ^24 ^^3 41 / dp e-'^^^-^-^) 

1=1 aibi 711=1 712=1 n3=0 n4=0 " 



(27r)5 



" dip 



X 



dk P Fnm{p,p;x) 



Df{K) + (1 - x') cos(/? Z^f (K) Z^f (ir) + (1 - x') cos ^Df{K) 



(124) 

The infrared singularity for the longitudinal component is cut off by the Debye mass, but 
that for the transverse component is regularized by the dynamical screening mass, ~ xiv?. 
The is infrared safe for a; 7^ 0, and is finite after integration over x. Consider the following 
integral: 



dk 



1, / 

In 



{k"^ + m\){k'^ + 4 \m\im? 



(125) 



2 g^N, ' 

where we use and ~ g'^N^T'^, and we keep only the logarithmic term. This expression 
is valid as long as trL 7^ 0. Under these approximations, the integral over k in Eq. (124) 
can be done, giving 

I,. .0, 1 



dk k^ \Dl + (1 - x^) cosif Dt\ ^ -(1 - cosy^)^ In 



g'N, 



(126) 



We can then perform the angular integral. 



dx 



" dip 
.2^ 



Fnm{P,P,x) (1 - cosy?)' 



XniP)Xmip) + 2 (p'^Xnip))' {p^Xm{p))' + {P^P') 



(127) 
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n+m 



P 



1=0 



+ 



where 



Cnm,l = E E ( 1 + (■S^ + 2/) j 6s+t,l 

s=0 i=0 ^ ^ 



128) 



This integral is symmetric under exchange of p and p', so we can replace + p'^ by 2p'^. 
Using 



dp p e 



2 „-/3p(ni+n3) 



48 r(5 + /) 



(rii +n3)3(n2 + n4)5+'r(5) 



(129) 
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20 4/in4 oooooooo 

r = ^ In 



i=\ ai,bi ni=l 7i2=l "3=0 7i4=0 H '?n'l 

n+m 



the matrix element £nm becomes 

\^^2]\rl n 

This matrix element is ~ T^g'^\\i{l/ g'^N^, so that rj ~ /{g'^\al/{g'^N^), as expected. 
Thus in the semi-QGP, the factors of temperature and the coupling constant are the same 
as in the perturbative QGP; the dependence upon the background field appears only through 
Q{Q,n). 

The next step is to sum over color. There are sixteen diagrams contributing to the 
viscosity. Interference terms are not planar diagrams, and so can be dropped at large A^^^. 
In Fig (2), the contribution of diagrams (a) and (6) coincides with that of {d) and (c). The 
color structure of diagrams (a) and (6) are 

Gia){Q.n) = exp (Q"(ni + ^3) + Q\n2 - rii) + Q'im - ^3) - Q'^{n2 + n^))] , 
e?(fe)(Q, n) = exp [t(3{Q''{ni + n^) + Q\n2 + n^) - g^(ni + rig) - Q''{n2 + n^))] , (131) 

respectively. Introducing 

mi = ni + ris , 1712 = n2 + n4, 1713 = 123, m4 = 114 , (132) 
the matrix element becomes 

_ lOT«(^^ivJX^ f 1 \ ^ ^ ^ ^""^im, n5+i) ^ 

TT^ WNj 2^ 2^ 2^ 2^ 2^ ml m^+' r(5) (...^ 

mi = l m2 = l m3=l m4=l ^=0 ^ 2 V / (13oj 

X (■^m4— m3-^m2+»n3— mi— m4 ~l~ -^mi— m3+m4-^m2— m4+m3) • 

This is the final form for Cnm- If the background vanishes, 

^3 C(2) C(4) In J ^ , (134) 

where ^(-2^) is the Riemann zeta-function. In particular, for the inverse of £, the zero- zero 
component dominates, 

= y ° C(2) C(4) lu (^) . (135) 
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FIG. 4: Matrix size dependence of Cqq[C ^]oo — 1- 



where C(2) = ttVG and C(4) = 7r790. At finite iV^, A^,^ is replaced by Nl{Nl - 1). In zero 



background field, 



' C(5) 



(136) 



where iV^ is replaced by (A^"^ — 1) at finite A^. 

In zero background field, at large A^ the shear viscosity is 

8T" 



"o 



45 Cooii 



540 C'(5) 



2\5 



(137) 



vr/ {gm,)nnl/{gm. 



This result was first obtained by Baym et al. [11]. Our result is valid at large Ac, with the 
correct result at finite is obtained by replacing the overall factor of A"^ in the numerator 
by A'2 - 1. 

Comparing with Eq. (54), the coefficient = 26.98. This result is very close to the exact 
result, = 27.13, which was obtained numerically by Arnold, Moore, and Yaffe [16]. We 
confirmed their results with several large matrices. In Fig (4) we show the dependence of 
[£^^]oo(^ = 0) upon the size of the matrix. The solution converges quickly as the size of the 
matrix increases, with the solution for one matrix element within ~ 0.6% of the exact result. 
Therefore, in the following we limit ourselves to the approximation of one matrix element. 

In order to compare the viscosity in the semi-QGP to report, we calculate the ratio TZ^i) 
in Eq. (53), 

nifl.^.'M^I^I^^. (138) 

Vpert («)]00 
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Inserting the explicit form, we obtain 



C(5) ^ 

^ ' n=l 



oo oo mi m2 

E\ ^ \ ^ <^mi 

Z-^ Z-^ Z-^ m'^ 



2C(4)C(2) ^ ^ ^ ^ 

-5V 7 ^2=1 m3=l m4=l ^ ^ 

1 -1 

^ (■^m4— m3^m2+m3— mi— m4 ~l~ -^mi — m3+m4^m2— m4+m3 ) 



(139) 



This is our main resuh. If the background field vanishes, i.e., all = 1, this ratio is unity, 
TZ^i) = 1. Conversely, if all in = 0, this ratio vanishes, IZi^i) = 0. Numerically we will see 
how TZ{i) changes as a function of i in the next section. Before going to the numerics, we 
estimate C near T^, where the Polyakov loop is small. Neglecting higher powers of i, and 
higher moments, in for n > 2, 

'^w = 270^(5) rrp-^-^^rr?- 

Thus the shear viscosity is ~ i"^ at small i. 

Let us see why the shear viscosity becomes small. Comparing Eqs. (31) and (35) to 
Eq. (116), one identifies 

^Uab ~ T^do{i) ; ^ nabncd(Tab;cd ~ ^^^f^ ■ (141) 

ab ab,cd 

As mentioned in Sec. Ill A, both the gluon number do ~ N^i"^, and the transport cross 
section multiplied by density squared. J2ab cd^abncd<7ab;cd ~ N^i"^, are small when i is. 

This suggests two important points. First, gluons (and other colored fields) are suppressed 
as T — )■ T^. This "bleaching" of colored fields is most natural in the confined phase; it is 
universal, and depends only upon the color representation, independent of the fields spin, 
flavor, or mass. Secondly, even if the coupling constant is moderate, the correlation between 
gluons in the semi-QGP effectively becomes large; this is why the viscosity is small when 
i is. Conversely, the interaction between small color representations becomes large: e.g., 
quark and antiquark scattering includes two channels, an adjoint an a singlet. For small i, 
the singlet channel dominates that in the adjoint. This reflects that as one approaches the 
confined phase, interactions are dominated by the formation of color singlet states, such as 
glueballs and mesons. 



29 




1.4 I ^ 

Gross- Witten 
1.2 - Step function 







0.2 



0.4 



0.6 



0.8 



1.0 



FIG. 5: Numerical results for TZ{i), Eq. (139). The two ansatzes for the eigenvalues are a step 
function (dashed line) and Gross- Witten (solid line). The dotted line denotes a simple analytic 
approximation at small i, Eq. (140). 

IV. NUMERICAL RESULTS 

A. The shear viscosity in the semi-QGP 

In this section we compute numerically the ratio TZ{i) in Eq. (139). This is the ratio of 
the shear viscosity in the semi-QGP, when (£„) ^ 1, to that in the perturbative QGP, when 
all (in) = 1- The results are plotted in Fig. (5) as a function of (i). Of course the ratio TZ{i) 
is a function not only of the first moment of the eigenvalue distribution, {£) , but of higher 
moments, The problem is that while {£) is easily extracted from numerical simulations 
on the lattice, the full eigenvalue distribution is not. 

To illustrate the range of possible results, as mentioned in Sec. II we consider two different 
distributions, which we hope are representative. The first is a step function, Eq. (22); the 
second is that of a Gross- Witten model, Eq. (24). The different results for 7l{i) is illustrated 
in Fig. (5): there the result for a step function is shown by a dashed green line, and that for 
a Gross-Witten model, a solid red line. The difference between the two distributions are at 
most a few percent, over the entire range of i. In fact, up to £ = 0.6, an approximate solution 
without higher moments in Eq. (140), denoted by a blue dotted line shown in Fig. (5), is 
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FIG. 6: Feynman diagram of quark contributions to the shear viscosity in the leading order. The 
straight hne denotes quarks or antiquarks, and the curly line denotes gluons. 

close to the other two distributions. This suggests that at least for the shear viscosity, the 
results are not very sensitive to higher moments of the eigenvalue distribution. We do not 
know if this is generic, or special to the shear viscosity, computed at leading logarithmic 
order 

We find suppression of the shear viscosity at small ~ £^ as £ — )■ 0. Conversely, we 

find that 7?.(£) is enhanced for I ~ 1.0. For £ ^ 0.9; 7?.(£) ~ 1.25; near £ = 1, an approximate 
form is ~ 1 + lATy/l — i. We do not have a simple explanation for this enhancement; if it 
is a peculiarity of working at leading logarithmic order, etc. 

We also computed the shear viscosity with dynamical quarks; details will be presented 
elsewhere [54]. As in the pure glue theory, the dominant processes to leading logarithmic 
order are those involving the exchange of a soft particle in the t-channel. The diagrams which 
contribute are illustrated in Fig. (6). There are contributions from the scattering between a 
quark and an antiquark, or scattering between two quarks. Fig. (6. a); scattering between a 
gluon and a quark. Fig. (6.b); and Compton scattering between quarks and gluons. Fig. (6.c) 
and Fig. (6.d). 

However, the physical processes which dominate are very different from the pure glue 
theory. At small i, the dominant scattering processes are not between two gluons, but 
between a quark and an antiquark, as shown in Fig. (6. a). Contributions to the shear 
viscosity arise from two parts, do and £oo; shown in Eq. (116). do is proportional to the 
number density; for gluons, this is of order i"^, while with quarks, it is ~ i. In £oo, gluon- 
gluon scattering gives £oo ~ the scattering of a quark and an antiquark includes a 
channel. Fig. (6. a), which is color singlet, so that £oo ~ Therefore, for a purely gluonic 
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FIG. 7: Flavor dependence of the viscosity. 

theory, the shear viscosity is r/ ~ dl/Coo ~ ~ with dynamical quarks, it is also 

rj ~ (Iq/Coo ~ {iy/i^ ~ i"^. The similar behavior of r] at small i must presumably reflect 
some more fundamental physics, although we do not know yet what it is. 

In the numerical results, we took the ratio of the number of flavors to color to be flxed at 
several values, equal to Nj/N^ = (the pure glue theory), 1/3, 2/3 and 1. The results are 
illustrated in Fig. (7) [30]. There is a weak, but non- negligible, dependence upon Nf/Nc. 

B. The ratio of the shear viscosity to the entropy 

In applications to hydrodynamics, what matters is not the shear viscosity by itself, but 
the ratio of the shear viscosity to the entropy density, rj/s. For A/" = 4 supersymmetric 
gauge theories, the gauge/gravity duality gives rj/s = l/iii [6, 7], which is conjectured to 
be a universal lower bound. 

In this subsection we compute the behavior of rj/s in the semi-QGP. In principle, with 
a complete theory of the semi-QGP the entropy density would follow from the eigenvalue 
density. Since such an eigenvalue density is not presently known, we use results directly 
from the lattice, siat [35]. This is not a signiflcant limitation, since any complete theory of 
the semi-QGP would necessarily give an entropy density in agreement with siat- 

We also need to know the temperature dependence of the coupling constant in QCD. We 
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FIG. 8: Ratio of the shear visosity to the lattice entropy. 

use the formula for the running of the couphng to two loop order, 

Ncg\T) 47r 



Nra, 



4- /3o(x)lnf + |Mlnlnf 

where the coefficients of the beta function at large are 

/3o=^(ll-2x), 
/3i=^(34-13x), 



(142) 



(143) 
(144) 



with X = Nf/Nc- We determine the renormalization mass scale, A, such that as{Tc) ~ 0.3 
[24]. This value of the coupling is moderate in strength, so that our computation, to leading 
logarithmic order, may not be reliable. To parametrize this uncertainty, we introduce a 
phenomenological parameter k to represent effects at higher order in the shear viscosity. 

In Fig. (8) we plot 77/siat, when = Nf, for three values of k, k = 8, 32, and 64. As 
the temperature decreases, the coupling increases, and the shear viscosity, r] ~ l/g'^{T), 
decreases. For larger values of k, the dependence on the temperature is weaker. 

The shear viscosity rj ~ /{g'^Nc)'^; with g'^Nc fixed at large Nc, rj ~ N^. Our compu- 
tation is vahd to leading logarithmic order, and includes the effects of quarks and gluons. 
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Hadrons, such as mesons and glueballs, are subleading, and contribute ~ 1 to 77. Thus when 
at small I ~ l/^c? our approximations break down. We expect that this is the reason why 
our results for r^/siat violates the conjectured lower bound, near and below T^. 

V. "BLEACHING" COLOR IN THE SEMI QUARK GLUON PLASMA 

We conclude with some general comments about our results, and some speculative re- 
marks. 

The simplest way of viewing deconfinement in a non-Abelian gauge theory is in analogy 
to ionization in an Abelian plasma [30, 45]. Define the lower limit of the semi-QGP as T~^^:^, 
and the upper limit, as T^^:^. Recall that in QCD, lattice simulations indicate that the semi- 
QGP exists between T~^-^ ^ 0.8 Tc and Tg^^^j ^ 3 — 4Tc, as discussed in the Introduction and 
Refs. [30-32]. 

There is no ionization of color below T~^^:^, so the only states are colorless hadrons. The 
ionization of color is partial in the semi-QGP, between T~j^j and Tg^^i' ^^"^ total above 
Tg^jjjj, in the perturbative QGP. A simple corollary of this picture is that scattering process 
involving colored states disappears as T — )■ T"^^; from above. One may say that the scattering 
of colored particles is "bleached" by the vanishing of the Polyakov loop. In this paper we 
computed how the shear viscosity is bleached, an example especially relevant for experiment 
[4]. 

One result which followed from our analysis, but which we did not stress, is the following. 
The Polyakov loop is related to the propagator of an infinitely heavy, test quark. Yet we find 
that the bleaching of dynamical, light colored fields, such as gluons and massless quarks, is 
identical to that of heavy quarks. That is, at least within the approximations in which we 
work, the bleaching of color is universal, independent of the mass of the dynamical fields. 

Thus it is necessary to discuss what assumptions are implicit in our analysis; these are 
equivalent to the statement that the semi-QGP exists only in a relatively narrow window 
about Tc. 

That Tg~jj^j is not too far below implies that it is reasonable to speak of confinement in 
a theory with dynamical quarks. It is not difficult to think of gauge theories for which this is 
not true. Consider a gauge theory with three colors and eight massless flavors, where lattice 
simulations appear to indicate that chiral symmetry is spontaneously broken in the vacuum 
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[55]. If so, then there is a chiral phase transition at which this symmetry is restored at a 
finite, nonzero temperature. Nevertheless, in this theory the chiral phase transition bears 
little relation to the deconfining phase transition of the pure gauge theory. This follows 
simply by counting the number of degrees of freedom: the hadronic pressure, from an ideal 
gas of 63 Goldstone bosons, is large, ^2/3 that of the ideal QGP. Concomitantly, presumably 
the expectation value of the (renormalized, fundamental representation) Polyakov loop is 
also significant at temperatures well below that for the restoration of chiral symmetry. In this 
case, one might say that the semi-QGP begins at a temperature far below Tc, but the physics 
is really just eight fiavors of quarks throughly wash out the Z{3) global symmetry of the 
pure glue SU{3) theory. The physics is dominated by chiral symmetry, and its restoration, 
with deconfinement a minor perturbation on that. 

The other assumption is that T^^^ is not much greater than Tc. Assume that the contrary 
were true, Tg^^i ^ '^c- then we expect that both the ratio of the pressure, to the ideal 
gas pressure, and also the (renormalized, fundamental representation) Polyakov loop would 
deviate from unity, from T^^^-^ all of the way down to Tc. While straight Polyakov loops 
presumably dominate near T^^^, it is unreasonable to expect that they dominate all of the 
way down to Tc, if T^^^ ^ Tc. Instead, we suggest that in this case it is necessary to include 
thermal Wilson lines, and the corresponding Polyakov loops, which oscillate an even number 
of times in r : — ?■ 1/T. (The number of oscillations must be even because a Wilson loop, as 
a quantity formed from bosonic fields, must be periodic in r.) As one goes to temperatures 
much lower than T^^^, the effects of oscillatory Polyakov loops become more significant. 
An effective theory of Polyakov loops can still be constructed; it is just that a new type of 
Polyakov loops have to be folded in. 

If effects from oscillatory Polyakov loops are important, then it is clear that the propa- 
gation of light fields would differ from that of heavy fields. A heavy field propagates in a 
straight line in imaginary time, but a light field performs a random walk; the lighter the 
field, the more the dominant paths in the path integral include those which fiuctuate from a 
straight path. Thus if T^^^ were much larger than Tc, the light fields would feel the effects 
of oscillatory Wilson lines, but the heavier fields would not. 

In fact, even in a pure gauge theory, where the semi-QGP exists only in a narrow region 
in temperature, there is no reason why oscillatory Polyakov loops could not be constructed 
and measured. On the lattice, it will be awkward to discretize them, but they could be 
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measured, on a sufficiently fine lattice. Our principal assumption is then that they are not 
needed when the semi-QGP exists in a narrow region in temperature. 

We comment that there is a soluble limit in which one can show that only straight 
Polyakov loops control deconfinement. For a theory on a sphere in the three spatial di- 
mensions, the radius R of the sphere can be made of femtometer dimensions, so that the 
non-Abelian coupling constant runs to small values. This theory is soluble, and has a true 
phase transition at infinite [56]. Near T^, the deconfining transition is controlled by the 
simplest Polyakov loop, a straight lop in the fundamental representation. At zero coupling, 
deconfinement occurs at a temperature T^, which is equal to a pure number times 1/R; at 
this point, the mass term for this Polyakov loop, alone, changes sign, jumping to a value of 
1/2. Although the full behavior of this model has not yet been computed, since the only 
dimensional scale in the problem is the radius of the sphere, it is most likely the semi-QGP 
only persists up to temperatures a few times l/R; i.e., T^^jTc is a number of order one. 

There is a heuristic explanation as to why lattice simulations appear to find that T^^^ 
is just a few times Tc. If T^^-^ ^ T^, then given how the couphng constant in QCD runs 
with temperature [24] , there would have to be a perturbative mechanism for generating the 
eigenvalue repulsion necessary in the semi-QGP [26]. This would indicate a perturbative 
instability for QCD in a thermal bath, which seems unnatural. Instead, lattice simulations 
find that the deviations from conformality, apart from the usual perturbative corrections 
from the conformal anomaly, are nonperturbative, due to corrections ~ [26, 34-36]. If 
so, in a narrow regime in temperature it is most natural that there are only a few operators 
which contribute, and that these few only involve straight Polyakov loops. 

Thus our conclusions about the universality of the bleaching of color are special to a 
semi-QGP which exists only in a narrow region of temperature. We also do not imply that 
all properties of fields in the semi-QGP are independent of their mass: only the bleaching 
of color. This result may be of significance for experiment. One of the real puzzles of the 
experimental data from the Relativistic Heavy Ion Collider is that the behavior of heavy 
quarks, such as charm, appears to be rather similar to that of light quarks: see Refs. [1, 2], 
and especially Ref. [3]. This is very difficult to understand if the behavior arises from 
energy loss, which is very different for light quarks than for heavy. Our analysis suggests 
a completely different mechanism may be responsible: the bleaching of color. Of course a 
more careful analysis is necessary in order to confirm this suggestion. 
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If RHIC is in a conformally invariant regime, as suggested by A/" = 4 supersymmetric 
gauge theories, then results for heavy ions at the LHC should be similar to RHIC: rj/s 
remains small, essentially unchanged with temperature, and the behavior of heavy quarks 
will remain like that of light quarks. 

In contrast, if the semi-QGP is valid, then there will be dramatic differences between 
RHIC and the LHC. If experiments at RHIC probe a region near Tc, then those at the LHC 
should probe temperatures significantly, perhaps a factor of two, higher. At the LHC, the 
ratio of the shear viscosity, to the entropy density, will be large. Also, in a perturbative 
QGP, the behavior of heavy quarks should differ significantly from that of light quarks. The 
difficulty is that these differences are only for initial times and temperatures; inevitably, a 
system at the LHC, even if it starts in the perturbative QGP, cools through the semi-QGP. 

In the end, we eagerly await the experimental results for heavy ions from the LHC, which 
will decide which theory is correct. 
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